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ABSTRACT 


The object of the thesis is development of a method for 
describing the behavior of a ship's material condition over 
time so that the way in which various maintenance strategies 
affect this material condition may be examined. The thesis 
demonstrates the manner in which these strategies may then be 
evaluated based on total expected cost of a fleet of ships. 


Ship failure and deterioration behavior is modeled as a 
transient Markov process, with a set of increasingly complex 
models. The most comprehensive of these allows such behavior 
to depend on the ship's present condition, its age, and the 
time since it was last overhauled. In the models a ship is 
treated as a single component, rather than an aggregation of 
several smaller components, and it is assumed that the ship's 
material condition can be categorized in a discrete manner. 


Use of the models is demonstrated by the use of two 
example cases, one of which is based on actual maintenance 
data for Navy ships. 


Thesis Supervisor: J. W. Devanney III 
Title: Associate Professor of Marine Systems 
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I. INTRODUCTION 


A. OBJECTIVES: 


It is the sm. of this thesis to develop an approach 
for examining the behavior of a ship's material condition over 
time and the manner in which that behavior is affected by var- 
ious alternative maintenance policies. We will do this under 
the premise that a trade-off exists between the investment 
cost of a fleet of ships and maintenance expenditures for the 
fleet under various maintenance policies, given an operational 
requirement of the form that some number of ships x must be 
available for operations with probability p. For a specified 
x and p, our goal will be to determine that combination of 
ship acquisition and overhaul policies with the lowest expec- 


ted cost, which obeys this constraint. 





B. BACKGROUND: 


In recent years nearly 20 percent, or approximately a 
billion dollars a year, of the Navy's entire operation and 
maintenance budget has been spent on the maintenance of the 


! This does not include the cost of naval man- 


Navy's ships. 
power, the cost of support ships and advance support bases, 
nor the capital cost of naval shipyards, all of which contri- 
bute significantly to the real cost of ship maintenance. Nor 
does it include any additional investment in ships that may 
result from choosing one maintenance policy over another. 

A significant part of the cost of ship maintenance is 
incurred during the ship's regular overhaul. A regular over- 
haul is a scheduled period (under present Navy policy), occur- 
ring at specified intervals, during which a ship is assigned 
to a shipyard for work of a relatively major nature. It 
includes maintenance, both corrective and preventive, as well 
as (usually) some items of modernization. 

Because a Significant portion of the maintenance dollar 
is spent on regularly scheduled ship overhauls, a great deal 
of attention has been focused on improving the effectiveness 
of overhauls, i.e., attempting to get maximum benefit from the 

lSee U.S. Office of Management and Budget, Budget of the 
United States Government and appendices thereto, 1969-1972, 
Government Printing Office, Washington, DC; and Department of 


the Navy, Naval Ship Systems Command, The Cost Estimate Fact 
Book, Volume I: Systems and Methods, Washington, DC, 1969. 
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work accomplished relative to the amount of time and money 
invested. This effort has generally proceeded along two 
lines: (1) By attempting to shorten overhauls while main- 
taining the same amount of work accomplished by improved 
planning, skill utilization, material procurement, etc., by 
such efforts as the Compressed Regular Overhaul (CROH) Pro- 
gram; (2) By examining the appropriateness of the length 
between regular overhauls and the manner in which a ship's 
material condition deteriorates between regular overhauls. 
It is in this latter area, perhaps, that the most effort has 
been exerted. 

In this area, again, two basic approaches seem to have 
been taken: The first of these is to examine separately the 
failure behavior of the various components which make up a 
ship and then attempt to synthesize this into a description of 
the behavior of the material condition of the overall ship. 

To conduct such an analysis for just one component, assuming 
the data to support the analysis is available, is a rather 
extensive task. Frequently, sufficiently complete or accurate 
data to adequately support such an analysis is not available 
for the component. In considering an entire ship made up of 
Several hundred subsystems the magnitude of the job becomes 
quite formidable. In the Navy, this is compounded by the lack 
of commonality between components of ships of the same type, 


or even of the same class of ship within a type. This has 
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resulted in part from the Navy's policy (in the past, at 
least) of letting contracts for ships of a given class to sev- 
eral different building yards without requiring commonality of 
equipment. 

A more promising path to take (and the second basic 
approach to the problem), at least in the short run, would 
seem to be to treat the whole ship as a single component. 
Indeed, using this approach, Farrar and Apple? were reason- 
ably successful in developing useful formulas for repair 
requirements as well as information on the difference between 
scheduled and unscheduled repairs. This was done in spite of 
much of their data being incomplete and somewhat inaccurate. 
Their work tends to support the basic feasibility of this 
second approach, in that a great deal of useful, yet relative- 
ly inexpensive information can be derived by considering the 


ship as a single unit. 


¿See Donald E. Farrar and Robert E. Apple, "Some Factors 
That Affect the Overhaul Cost of Ships: An Exercise in Sta- 
tistical Cost Analysis," Naval Research Logistics Quarterly, 
vol. 10, no. 4, December 1963, pp. 335-368, and "Economic Con- 
siderations in Establishing an Overhaul Cycle for Ships: An 
Empirical Analysis," Institute of Naval Studies Research Con- 
tribution No. 7, April 1964, Naval Engineer's Journal, vol. 
77, no. 1, February 1965, pp. 69-78. 
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C. THE APPROACH: 


Here we pursue the "whole ship" approach, and develop an 
increasingly general set of stochastic models describing the 
behavior of a ship's overall material condition as the ship 
ages and as alternate overhaul strategies are employed. These 
models will be designed to be intuitively appealing to, and 
in fact to make use of, the informed judgment of Navy decision 
makers. A model which a decision maker does not understand, 
or which he finds at odds with his intuition, is unlikely to 
see much application. 

Once a model has been developed through which a reason- 
able approximation of the behavior of a ship's material con- 
dition over its life can be made, a number of different 
investment and maintenance policies will be studied for their 
effect on the total cost of meeting a specified operational 
requirement. In particular, once the behavior of a single 
ship has been modeled, information can then be determined for 
an entire fleet of such ships. For example, the expected num- 
ber of ships in overhaul at a given time for a given fleet 
Size or, given the requirement for maintaining a specified 
number of operational ships with a specified probability, the 
required fleet size can be determined. We will show this in 
the thesis. 

In this thesis we will consider cost as belonging to two 


categories: overhaul cost, and the ship acquisition and life 
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cycle costs not including overhaul cost (reduced to present 
value). The latter will henceforth be referred to as "ship 
cost". In order бо be valid for decision making, this must 
include the present value of the ship's entire life cycle 
costs less overhaul costs, including cost of crew, etc., as 
well as the actual acquisition cost. If we know or can reas- 
onably estimate ship cost as well as overhaul cost, and if we 
can determine the expected number of overhauls the ship will 
experience under a given policy, then we can determine total 
expected lifetime cost. 

The models developed in the thesis will include the capa- 
bility to investigate the effect on total fleet cost of var- 
ious overhaul policies and the required fleet size (number of 
ships) resulting from these policies. 

We will model ship failure and deterioration behavior as 
a transient Markov process. Beginning with an extremely sim- 
ple process we will elaborate it until it has the ability to 
reflect observed behavior in a reasonably complete manner. 
Once this has been done, data extracted from actual ship main- 
tenance records will be introduced and use of the most com- 
prehensive model as an aid to decision making will be demon- 
сва сеа. 

This approach to modeling a ship's change in material 
condition over time represents, it is felt, a somewhat unique 


conceptualization of the problem. Although Farrar and Apple 
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treated the ship as a single component, their approach depen- 
ded primarily on the use of regression analysis. 

The model should be applicable to both naval and commer- 
cial ships, at least in basic concept. However, in the 
thesis, where specifics are required in either developing or 
demonstrating the model, they will be derived from and geared 


to destroyer-type ships of the U.S. Navy. 
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D. SUMMARY: 


The remaining chapters of the thesis are summarized as 
follows: In Chapter II we first establish the underlying 
assumptions for the model. We then develop a simple model in 
which two categories of ship material condition are assumed 
(consider these as "failed" and "not failed") and in which 
the probability of failure is dependent solely on the ship's 
age. Chapter II concludes with a hypothetical example demon- 
strating use of the model as developed at that point. 

In Chapter III we introduce into the model the assumption 
that not only will the probability of failure change with 
ship age, but with the length of time since the ship was last 
overhauled as well. We then expand our binary failure assump- 
tion (two categories of material condition) to one which 
allows for four categories of ship material condition, in 
which the worst of these, condition 4, corresponds to failure 
in the binary model. We conclude Chapter III with an example, 
using our four condition model, based on data extracted from 
the maintenance records of the Commander, Cruiser-Destroyer 
Force, Atlantic (COMCRUDESLANT). 

Chapters II and III of necessity deal to a certain extent 
with the mechanics of computation for the model. It is not 
until we reach Chapter IV that we are actually able to apply 
the model to a comparison of alternative overhaul policies, 


utilizing the example of Chapter III as a basis. In Chapter 
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IV we examine the effect of overhaul policy on the number of 
ships necessary to meet a given operational requirement with 
a specified probability. Chapter IV compares overhaul poli- 
cies of several specified intervals with the policy of over- 
hauling a ship only upon failure, i.e. a "demand" dependent 
policy. Chapter IV further considers the effect of varying 
the specified operational requirement and the effect of im- 
puting some penalty cost to an unscheduled overhaul, i.e., 
one which occurs on demand due to a ship failure. 

In Chapter V, after discussing some possible directions 
for further work in this area, we arrive at the following 


conclusions resulting from the thesis: 


(1) That the behavior of a ship's material condition 
over time can be successfully modeled by a transient Markov 


process. 


(2) That such a model can be of use to Navy decision 


makers in comparing alternative overhaul policies. 


(3) That actual implementation of the model will require 


initially a major data collection and reduction effort. 
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II. DEVELOPMENT OF THE BASIC MODEL 
А. ASSUMPTIONS: 


The basic assumptions underlying the models used in this 
report are: (1) that the change in a ship's material condi- 
tion over time may be described by a Markov process, (2) that 
the ship can be classified into one of two or more discrete 
categories of ship material condition, and (3) that some 
reasonable approximation of the probabilities of remaining in 
the same or changing to a different material condition cate- 
gory as the ship ages can be determined. 

A first order Markov process is one in which only ве 
last state occupied by the process is relevant in determining 
its future behavior. Under the Markovian assumption the prob- 
ability of making a transition to each state of the process 
depends only on the state presently occupied.? 

As noted earlier, a ship, especially a warship, is suf- 
ficiently complex that its overall material condition is com- 
prised of the net effect of the material condition of a large 
number of component systems and subsystems. Further, it can 
be significantly affected by a number of diverse factors, such 
as the intensity of operations, the funds available for re- 


pairs and spare parts, the qualification and motivation of the 


’See Ronald A. Howard, Dynamic Probabilistic Systems, Vol- 


ume I: Markov Models, John Wiley and Sons, New York, 1971,p.3. 
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operating crew, etc. Consequently, a means of precisely 
determining a ship's overall material condition has thus far 
been somewhat elusive. Efforts to devise some sort of 
material condition index have been initiated and are contin- 
uing; however, no such usable index exists at this time. 

It seems rather that the best means to try to determine 
a ship's overall material condition is to ask the operator/ 
maintainer, in the present case the type commander. Such an 
approach is necessarily subjective and is more likely to be 
qualitative than quantitative in nature. The type commander, 
supported by his staff, is nonetheless likely to have a good 
"feel" for the condition of his ships since he is in frequent 
contact with them and has the benefit of experience with var- 
ious ships. Based on this, the model assumes that a type com- 
mander, utilizing the collective knowledge and experience of 
his operating and maintenance staffs, can (and in practice 
does, although perhaps not formally) categorize his ships into 
two, three or four different levels of material condition. 
Some artificiality exists in that the model assumes discrete- 
ness where in fact a continuum exists. Nonetheless, exper- 
ience has indicated that a type commander can fairly well tell 
which of his ships he considers "good", which are "bad", and 
which he would group in the middle as not belonging to either 
extreme category. 

Initial estimates of the probability of transition from 


one of these categories to another, as the ship ages, must be 
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made. Obviously the best means of doing this is to observe 
the behavior of a large number of similar ships over a period 
of time. Unfortunately, in the Navy to date there has been 
little or no formal categorization of total ship material con- 
dition as described above, so no past records are available 
from which estimators of transition probabilities may be 
determined directly, although a similar system is in effect in 
the Navy to denote the seriousness to the mos Operating 
capability of material casualties to individual components 
(part of the Casualty Reporting, or CASREPT, system).  Conse- 
quently, the best approach would seem to be Bayesian in 
nature, where prior estimates of the behavior of a ship's 
material condition are established based on the best data and 
informed opinion available, and then are updated as better 
information becomes available. A discussion of how such prior 


estimates might be made will appear later. 
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B. THE TWO CONDITION, AGE DEPENDENT MODEL: 


Given the foregoing assumptions, we will begin our ana- 
lysis by considering a relatively simple case. In this sim- 
plest case, two possible conditions are assumed: these can 
be considered "good" and "bad", "operating" and "failed", etc. 
They will be designated as "1" and "2", respectively. 

Intuitively, it seems reasonable to expect a higher prob- 
ability of failure for an old ship than for a new one. While 
this may or may not be true, it would be desirable if the 
model were able to accomodate such a possibility. In order 
adequately to describe a state, therefore, both the ship's 
condition and age must be specified. Let the state indices be 
age t and condition k; the state then is described by the 


ordered pair: 


(k,t) = ship is t units old and in condition k, 
КЕК, ТЕН, 
where K = {1,2} for a two condition model, 


andern — про а с T being the ea ËN 
age the ship will reach) of the ship in 
terms of the basic time unit. 

The assumption is made that the ship will occupy the 
state for all of a given time increment. This should be reas- 
onably valid if the time increment is relatively short. In 
fact, it could correspond to a periodic appraisal of the con- 


dition of a fleet of ships, so that transition would take 
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place only when the ships' condition was reappraised. In this 
simple model, it is assumed that the decision maker, upon 
finding a ship in condition 2, would immediately take correc- 
tive action consisting of an "overhaul". In other words, а 
"demand" (for overhaul) condition would be said to exist any 
time the ship reached condition 2. The objective of over- 
hauling a ship presumably is to improve its material condi- 
tion. It will therefore be assumed that the overhaul will 
with certainty return the ship to condition 1. 227 2:27” 

Since an overhaul, under the latest assumption, in effect 
represents a transition from condition 2 to condition 1l, it is 
convenient to take as the basic time unit a period equal to 
the approximate length of overhaul of a ship, or a period of 
which length of overhaul is a multiple. A review of 59 recent 
overhauls of a variety of COMCRUDESLANT destroyer-type ships 
indicated an approximate four month overhaul length, on the 
average. The average time between overhauls was approximately 
40 months. For this model therefore, unless otherwise stated, 
the basic unit of time is assumed to be four months, although 
the actual length of the increment should in general be of 
little importance to the applicability of the model. Four 
months should, however, provide sufficient time for the deci- 
sion maker to perceive a change in condition. It provides, 


for a ship life of, say, 20 years, some 60 different ages. 


ЧЕ sn 1790 


и , 
7 


A 
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Transition probabilities will be described by p(kli,t), 
which is the оао of the ship's being in condition k at 
age t+1 given that it was in condition i at age t. What 
results, then,is a transient Markov process lasting the dura- 
tion of the ship's life (T time units) in which a state has 
zero probability of occupancy more than once. 

With this model, multi-step transitional probabilities 
may be calculated. If some initial state is fixed, the multi- 
step transition probabilities will provide the probability of 
entering a given state (i.e., a given condition at a given 
age). Multi-step transition probabilities will be described 
by P(k,t|k б.) which is interpreted as Príentering state 
(k,t)|process began in state (k ,t.)). Since at each transi- 
tion t must increase by exactly 1, the number of steps in the 
transition is given simply by EE For this model, multi- 
step transition probabilities are calculated by the recursion 


formula: 


2 
P(k,tik st) Sal Ровер око) РАМ код) 
k = 1,2. 


Of particular interest in this model will be the probability 
of entering a given state, given that the ship is new as an 
initial condition. If "new" is assumed to mean that the ship 
is in condition 1 at time zero (in other words that the ini- 


tial state is (1,0)), then this probability would be described 
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by P(k,tl1,0) for the probability of entering state (k,t). 
For ease of notation henceforth this will be written simply 
as P(k,t). 

This model will provide multistep transition probabili- 
ties (as has been described) and therefore the probability of 
the ship's being in overhaul at any given age. From this the 
expected number of overhauls over the ship's life may be 
determined. This is simply the probability of being in over- 
haul in each time unit of a ship's life, summed over the 
ship's life. Since we are assuming that a ship found to be in 
condition 2 will automatically undergo overhaul, we may write: 
E(number of overhauls over ship's life) - I0. The 
upper limit of T-l instead of T reflects e assumption that 
no ship would be overhauled during the last time unit of its 
life, just in time to be removed from service. 

If we desire to determine the probability distribution 
of the number of overhauls the ship will experience during 
its life (or up to a given point in its life), we must do the 
following: Define mT(k,t,x) as the probability of the ship's 
being in state (k,t) and having experienced x overhauls 
(including the current one, if k=2) given that the ship was 


new at t=0. This probability may also be calculated by means 


of a recursion formula: 
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т(1,6-1,хур(к | 1, #2) кеш, 
i 


И МО 


1 
mike) — 


T(i,t-1,x-1)p(k|i,t-1),k=2. 
1 


| 2 116) 


i 
We may then find the probability distribution for the number 
of overhauls in a ship's lifetime as follows: 


Pr{ship will experience 


2 
x overhauls in its lifetime} толк, т-а 


where T is the ship's life in terms of the basic time unit. 
We again use the value for t=T-1 instead of T. 

Now consider a fleet of identical ships, all of the same 
age. This is obviously an oversimplification, but will serve 
to demonstrate the manner in which the model may be used in 
considering groups of ships. Given the above assumptions and 
independence between failures of different ships, the proba- 
bility of having a certain number of ships in overhaul at any 
given time out of a fleet of a given size follows a binomial 
probability mass function. ТЕ а denotes the probability of 
a single ship being in overhaul at a given age, the proba- 
bility of having y ships of that age in overhaul out of a 


fleet size of N ships is given by: 
EES 


This is also, of course, the probability of having N-y ships 
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available for operations out of a fleet of N ships (assuming 
that a ship not in overhaul is available for operations). 

From this probability mass function can be determined the 
fleet size required to provide some number of operational 
ships (as a minimum) with a specified probability. Given such 
a requirement, a fleet size necessary to support it is as- 
sumed. The probability of having at least the minimum number 
of operational ships is computed using the binomial probabi- 
lity mass function. If this probability is greater than the 
specified probability, then the assumed fleet size is de- 
creased until the requirement is just met. If this proba- 
bility is less than that specified, the assumed fleet size is 
increased until the requirement is met. 

Since the expected number of overhauls over a ship's life 
may be determined from the model, assume that both ship cost 
and overhaul cost are known. Then the expected total cost 
may be determined for a single ship by multiplying the over- 
haul cost by the expected number of overhauls and adding the 
result to the ship cost. Once this cost has been calculated 
for a single ship it is known for a fleet (of a given size) of 


such ships. 
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С. АМ ЕХАМРЬЕ: 


Before proceeding with development of more comprehensive 
models it may be helpful to demonstrate application of the 
relatively simple model developed thus far in a hypothesized 
example. 

Suppose, for the example, that the life of the sample 
ship is expected to be 10 years, or 30 time units; further, 
suppose that the probability of "failure" increases in stages 


as the ship ages, as follows: 


ЕТ EE 
Щи: ое 
РИ Е) =. fae oe Ey 
.4 Е= E 
KÉ t-24,...,29 


It is assumed that the ship begins in state (1,0), and that 
if it should "fail", overhaul of the ship would begin 
immediately and would, with certainty, return the ship to 
condition l at the beginning of the next time unit. 

The following have been computed for the example: Values 
НРК, for all k and t; values of P(k,t|k,,t,) for t 
through ten time units; values of m(k,t,x) for all possible 
k, t and x; and the distribution and expected number of over- 
hauls the ship will experience in its life. All possible 


P(k,t|k_,t,) were not computed as the result becomes somewhat 
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unwieldy. For example, computing all values through 30 time 
units results in approximately nine times the output as for 
ten units. Because of this, and having demonstrated the 
manner in which multi-step transition probabilities may be 
calculated, henceforth only those values specifically 
required for some application will be computed, except for 


the special case of P(k,t) which is in fact the family of 


state probabilities, given that the ship began in the "new 
condition. 

The computer program used for the calculations was 
EEwoloped from the recursion relationships of the model. The 
program, together with the results for the example, comprise 
Appendix A to the thesis. The program is written in the PL/1 
programming language but the logic should present no problem 
to FORTRAN programmers. 

The expected number of overhauls the sample ship will 
experience over its life is 6.38. The probability distribu- 
tion of number of overhauls is shown in Table l. These prob- 
abilities are determined by summing 1(1,29,x) and 7(2,29,x) 
for all values of x. A ship age of 29 instead of 30 is used 
since it is assumed that it would be a waste of resources to 
overhaul a ship during tue final time unit of its life. The 
special family of transition probabilities P(k,t) are given 
by Table 2. Since these give the probability of being in a 
Specific state, given that the ship started "new", they will 


henceforth be referred to simply as "state probabilities" to 
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Number of Overhauls 


in Ship's Life Probability 
0 .00002 
1 . 00053 
2 .00499 
3 .02287 
4 ‚08295 
5 .17352 
6 .24346 
7 23183 
8 · 14980 
9 . 06494 
10 .01844 
11 .00328 
12 .00034 
13 „00001 
14 .00000+ 
15 .00000+ 


Note: 

Because of the certainty of returning to condition 1 
after overhaul, the ship cannot experience "failure" (and 
therefore overhaul) during two consecutive time units. There- 
fore the probability of more than 15 overhauls is zero. 


TABLE 1 
Probability Distribution of Number of 


Overhauls for the Two Condition, 
Age Dependent Model Example 





k Е P (b) 
1 1 . 900 
2 1 .100 
1 2 . 910 
2 2 .090 
1 3 . 91919 
2 3 .091 
1 4 . 909 
2 4 .091 
1 5 „909 
2 5 .091 
1 6 . 909 
2 6 .091 
1 7 . 818 
2 7 . 082 
1 8 . 836 
2 8 .164 
1 9 .833 
2 9 ‚167 
1 10 . 833 
2 10 .167 
1 11 „833 
2 11 ler 
1 12 . 833 
2 12 ‚167 
1 13 „750 
2 T3 .250 
1 14 лав 
2 14 25 
TABLE 2 


State Probabilities for the Two 
Condition, Age Dependent 
Model Example 
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N F џин NY мн мю ~ ~ [LS ~ KF NO FP мю PFP N F ~ ҥ N PFP ~ P N ka N P N F N ki 


t P(k,t) 
15 «767 
15 239 
16 . 770 
16 .230 
I7 „169 
1% 223] 
18 .769 
18 a! 
19 . 692 
19 308 
20 2723 
20 nove 
21 e 
21 - 299 
22 . 716 
22 .284 
28 .714 
23 .286 
24 .715 
24 . 285 
25 .643 
25 2:35 7 
26 . 679 
26 "921 
27 . 661 
27 .339 
28 . 670 
28 . 330 
29 . 665 
29 . 395 
30 . 667 
30 1999 


TABLE 2 (continued) 


30 
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differentiate them from other multi-step transition proba- 
bilities. 

Under the assumption that a ship found to be in condition 
2 will automatically — overhaul, the state probability 
P(2,t) is also the probability of the ship's being in overhaul 
at age t. For the present example, assume that the decision 
maker desires to maintain a fleet of ships identical to, and 
of the same age as, the sample ship. Assume further that he 
wishes to have a probability (say .95) of having a certain 
number of ships operational (that is, not in overhaul), say 
30 ships. The problem is to determine the size of fleet 
required to give a probability of at least .95 of 30 or more 
operational ships. 

With these assumptions, the required fleet size may be 
determined as follows: First, determine the highest probabil- 
ity of being in overhaul (that is, the maximum value of P(2,2) 
at any point in the ship's life. For the example this is 
P(2,25) (from Table 2) which has a value of .357. This repre- 
sents the worst possible situation. Next, using the binomial 
probability mass function, obtain the probability distribution 
for the number of operational ships using the value .357 and 
an assumed fleet size. Increase or decrease the fleet size as 
required until the requirement of a .95 probability of 30 
operational ships is just met. This is, of course, the smal- 


lest fleet size that will meet the specified requirement. 
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The minimum fleet size satisfactory to meet the require- 
ments specified in the example is 56. The computer program 
developed to accomplish this is contained in Appendix B. The 
initial fleet size assumed was 50 ships. The probability den- 
sity for the number of operational ships for a fleet size of 
56 is provided in Table 3. As would be expected, the proba- 
bility of 30 or more operational ships is .963 (slightly 
greater than .95 due to discreteness). 

It should be emphasized that use of the binomial proba- 
pility mass function is permitted only by the fact that the 
ships are identical and of the same age. The problem becomes 
more complex (in computation, although not in conceptualiza- 
tion) for a mixed fleet. 

Obviously, with the expected number of overhauls deter- 
mined, if a cost per overhaul can be estimated or is known, 
then overhaul costs over a ship's life can be а 
Assume the following representative costs for the example: 
Ship cost, $40 million; overhaul cost, $1.5 million. Then 


the mean of the total cost for the example is: 


$40.0 million + 6.38 x $1.5 million 


= $49.57 million. 


This can be expanded to include the entire fleet. 
Although the data is hypothetical and this model is rela- 
tively simple, the preceding example should serve to indicate 


the basic mechanics and some potential uses for such models. 





Number of 
Operational Ships Probability 
0 „ 0009004 
il .00000+ 
2 . 00000+ 
3 .00000+ ` 
4 .00000+ 
5 -00000+ 
6 .00000+ 
7 .00000+ 
8 .00000+ 
9 .00000+ 
10 . 00000+ 
11 . 00000+ 
12 .00000+ 
13 . 00000+ 
14 . 00000+ 
15 .00000+ 
16 . 00000+ 
17 .00000+ 
18 .00000+ 
19 .00000+ 
20 . 00001 
21 . 00003 
22 . 00008 
23 . 00022 
24 ‚00053 
25 . 00123 
26 . 00264 


TABLE 3 


Probability Distribution of the Number 
of Operational Ships for a Fleet of 
56 Ships (Two Condition, Age 
Dependent Example) 
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Number of 


Operational Ships 


2] 
28 
29 
30 
31 
32 
33 
34 
35 
36 
37 
38 
39 
40 
41 
42 
43 
44 
45 
46 
47 
48 
49 
50 
51 
52 
53 
54 
55 
56 


TABLE 3 (continued) 


Probabilit 


. 00527 
‚00982 
‚01706 
. 02762 
. 04167 
. 05857 
. 07664 


. 09327 


.10548 
.11070 


.10766 
.09685 


.08042 


. 06149 


. 04317 
.02774 
.01625 


. 00864 


. 00414 


. 00178 


. 00068 
. 00023 
.00007 


. 00002 


. 00000+ 
. 00000+ 


. 00000+ 


. 00000+ 


. 00000+ 
. 00000+ 
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The following chapter will deal with development of more com- 


prehensive models based on this one. 
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III. DEVELOPMENT OF MORE COMPREHENSIVE MODELS 
A. THE EFFECT OF OVERHAUL ON TRANSITION PROBABILITIES: 


As a first step in developing a more comprehensive model, 
consider the following: If one is willing to agree that some 
policy of overhauling ships on a regularly scheduled basis may 
be the most economical, one is implicity assuming that an 
overhaul may in some way "improve" succeeding transition prob- 
abilities. Otherwise, if no such improvement can be attained 
by an overhaul, there is no reason for overhauling. If a ship 
is just as likely to "fail" immediately following an overhaul, 
at some given age, as it is without a recent overhaul, at the 
same age, then there is obviously no justification for expen- 
ding the resources to overhaul the ship until it actually 
fails. 

It therefore appears that to describe adequately states 
and transition probabilities requires knowing the time (i.e., 
number of time units) since the last overhaul, as well as the 
ship's age and condition. Therefore, let the state indices 
now be age t, condition k, and time since last overhaul 3; 


the state is described by the ordered triplet: 


(k,t,j) = ship is t units old and has gone j units 
since its last overhaul, and is in con- 
dition k, keK, teH, jeJ, 


where again, K = {1,2} for a two condition model, 
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(1,2,...T), T being the life of the 21 


ra 
|| 


in terms of the basic time unit, 

and J = {1,2,...m}, m being the longest period 
a ship could possibly go without an 
overhaul; for a "demand" dependent 
(i.e., overhaul only upon ship's reach- 


ing condition 2) overhaul policy, m-T. 


For a ship that has yet to undergo its first overhaul, j will 
be taken as equal to t. 

Transition probabilities (one-step) will be described by 
p(k|i,t,j), which is the probability of the ship's being in 
condition k at age t+l, j+l units since its last overhaul, 
given that it was in condition i at age t, j units since 
last overhaul. 

It will now be necessary to describe multi-step transi- 


tion probabilities as follows: 


P(k,t,3|ko, t,:3,); 
Prlentering state (k,t,j) |transition began in 


State Коко“ Јо) Ј- 


Following the convention adopted earlier, the special family 
of transition probabilites Е ШОО О), the "state" proba- 
bilities, will be designated by P(k,t,j). 

The expression m(k,t,j,x) will be defined as the proba- 
bility of the ship's being in state (k,t,j) having experienced 


x overhauls.  Recursion formulas may be determined for multi- 
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step transition probabilities and values of m(k,t,j,x), given 
that the process started in state (1,0,0), for this two-condi- 
tion, age and time-since-overhaul dependent model. These 

have been developed in a manner parallelling those developed 


earlier and are shown in Appendix C. 
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B. THE MULTI-CONDITION, AGE DEPENDENT MODEL: 


While a structure which provides for only two categories 
of material condition serves sufficiently well to demonstrate 
Ae mechanics of the model, to say that a ship is either in 
good condition or is in bad condition and therefore in need 
of overhaul, with no middle ground, is intuitively unappeal- 
ing. ТЕ seems more likely that the decision maker could group 
his ships into three, four, or even more approximately dis- 
crete groups based on his subjective (but informed) opinion. 

At this point it will be assumed that a four condition 
model is reasonable in that it provides a wider range for 
categorizing a ship's condition without providing so many 
choices as to make such categorization overly difficult and 
arbitrary. There will be borderline cases, to be sure, re- 
gardless of the number of categories chosen. 

First, we will consider a four condition model in which 
transition probabilities are unaffected by the time since 
the ship was last overhauled (although they still can vary 
with ship age). In this model, the beneficial effect of an 
overhaul can be represented by a change to a better condition, 
e.g., condition 1 with certainty or condition 1 and condition 
2 with probabilities summing to unity, etc. Transition prob- 
abilities can be used to describe a gradually deteriorating 


material condition; all that is necessary is to structure the 
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transition probabilities so that, over time, the ship's con- 
dition is more likely to get worse than it is to get better. 
With this model, calculation of the entire set of (multi- 
step) transition probabilities remains reasonable, although 
still cumbersome; further, we can still determine, within 
reasonable computational limits, the probability distribution 
for the number of overhauls which the ship will experience 
in its life, as well as the expected value of this number. 
Computation of these probabilities is accomplished in 
a manner identical to that described in Chapter II, except 
that the set of possible conditions which the ship might 
teach, K, is {1,2,3,4}. 
Such a model will permit us to make comparisons between 
a demand dependent overhaul policy and one in which the ship 
is overhauled at regularly scheduled intervals (i.e., at cer- 
tain regularly spaced points in its life). Such a scheduled 
overhaul policy will have the advantage of permitting at least 
part of the overhauls in a ship's life to be scheduled in 
advance rather than resulting from the ship' reaching a demand 
condition (i.e., unscheduled). However, with this particular 
model, we cannot, of course, base a ship's entry into a 


scheduled overhaul on the time since it was last overhauled. 
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C. THE MULTI-CONDITION, AGE AND TIME-SINCE-OVERHAUL 


DEPENDENT MODEL: 


The most comprehensive model we will consider combines 
both the dependence on ship age and time since overhaul of 
the model of section A with the multi-condition capability of 
section B. This allows ship failure behavior to depend in an 
arbitrary manner on ship age, time since overhaul, and present 
material condition. This gain in generality is not without 
its sacrifice. As we will see, by requiring three indices 
to describe a state we have lost the ability to obtain some 
information about the process due to computational limita- 
tions. The recursion formulas for this multi-condition, ship 
age and time-since-overhaul dependent model are shown in 


Appendix C. 
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D. AN EXAMPLE FOR THE MOST COMPREHENSIVE MODEL: 


It follows from the description of this model that there 
are T x m sets of transition probabilities, since probabil- 
ity of transition is dependent on both ship age and time since 
last overhaul. If one considers a ship's life, T, to be about 
25 years, in the case where a ship isn't overhauled until the 
"demand" condition is reached, m is also 25 years. This rep- 
resents 75 time units in the model, which implies 75 x 75 = 
5625 possible sets of transition probabilities. This is 
clearly unmanageable. The answer, therefore, is to attempt 
to group ships which seem to follow similar "failure" rates 
into categories based upon age and time since overhaul. 

Following this approach, use of the multi-condition, 
age and time-since-overhaul dependent model will be demon- 
strated utilizing transition probabilites estimated from the 
RAV (restricted availability) history of 33 CRUDESLANT des- 
troyer-type ships. The manner in which these probabilities 
were estimated is described in Appendix D. The ships were 
grouped into three categories by age: Category I, 1 to 3 
years; category II, 4 to 13 years, and category III, 14 to 25 
years, inclusive. Three categories of time since overhaul 
were decided on: one (four month) time unit, two to three 
time units, and four or more time units. This results in nine 
sets of transition probability matrices. These are shown in 


Table 4. 
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Age Category I 


j = 1 
PN 205301505 
E 253 
0)3|1,t,1) = 0.121 
p(4]1,t,1) = 0.121 
p(i|k,t,1) are of no concern for i = 2,3,4, since ship can 
only be in condition 1 at 3 = 1. 
j = 2,3 
E t,J) 0.663 p(1 20000-900263 и = 0-52 
el2]1,t,j) 0.263 p(2|2,t,j) = 0.663 p(2| 3,67 10) =o, 263 
ESI, t,j) 0.037 p(3|2,t,j) = 0.037 p(3|3,t,j) = 0.437 
А |1,Е,3) 0.037 p(4|2,t,j) = 0.037 p(4|3,t,j) = 0.037 
J = 4 or greater 
ВИТ, Е, -) 0.682 212 = 0.282 p(1|3,t,)) 02235 
Er, t,j) 0.282 p(2|2,t,j) = 0.682 p(2|3,t,j) = 0.282 
ша, +, 3) 0.018 В =о ош p(3|3,t,j) = 0.418 
p(4|1,t,3j) 0.018 p(4|2,t,j) = 0.018 p(4|3,t,j) = 0.018 
TABLE 4 


One-Step Transition Probabilities for the 


Multi-Condition, Age and Time-Since- 
Overhaul Dependent Example 





Em t,j) 
p(2|1,t,]3) 
HSILt,j) 
(4 1 ,Е,3) 


E, t, 3) 
p(2[1,t,j) 
AS 1,t,)) 
p(4|1,t,3) 


0.662 
0.262 
0.038 
0.038 


0.594 
0.194 
0,106 
0.106 


Age Category TT 


j = 1 
ОЕК 30333 
EIFE OE y 
p(3|1,t,1) = 0.250 
p(4 1,t,1) = 0.250 

j = 2,3 

p(1|2,t,j) = 0.262 
BL 2) 2 ON 
p(3|2,t,j) = 0.038 
p(4|2,t,j) = 0.038 
1 + 4 or greater 
р(1|2,+,3) = 0.194 
p(2|2,t,j) = 0.594 
p(3|2,t,j) = 0.106 
p(4|2,t,j) = 0.106 
TABLE 4 (continued) 


ры Е 
p(2|3,t,3) 
p(3|3,t,3) 
Pie tea 


БОВЕ) 
P ZIE EJ) 
PO Е) 
p(4|3,t,j) 
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0202 
0.262 
0.438 
0.038 


0.194 
0.194 
0.506 
0.106 





Olti) 
E 1,t,j) 
p(3|1,t,)j) 
p(4|1,t,3) 


it) 
p(2|1,t,3) 
F |l,t,j) 
jt) 


0.664 
0.264 
0.036 
0.036 


0.627 
0.227 
0.073 
0.073 


Age Category TITI 
=н 


Ol | eens) 
p 21; Ui) 
p(3|1,t,1) 
p(4|1,t,1) 


0.571 
0.285 
0.072 
0.250 


ТАВІЕ 4 


j = 2,3 


p(1|2,t,3) 
p(2|2,t,j) 
p(3|2,t,j) 
р (4|2,+,9) 


p(1|2,t,3) 
p(2|2,t,3) 
p(3|2,t,3) 
p(4|2,t,j) 


0.264 
0.664 
0.036 
0.036 


3 = 4 or greater 


0.227 
0.627 
0.073 
0.073 


(continued) 


р, 
p(2[37 6/59 
aa 
pM, t. up 


p(1|3,t,3) 
p 225659) 
p(3|3,t,j) 
p(4|3,t,3) 
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0.264 
0.264 
0.436 
0.036 


0.227 
0.227 
0.473 
0.073 
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The computer model used for the calculations is shown in 
Appendix E. This model omits two features of the model used 
in earlier calculations (Appendix A). These are: (1) The 
computation of multi-step transition probabilities, for 
reasons discussed earlier, and (2) The computation of 
T(k,t,j,x), i.e., the probability of being in a given state 
and having experienced a given number of overhauls. The com- 
putation of n(k,t,j,x) vas omitted due to the computer space 
required to accomodate this four dimensional array. For the 
demand dependent policy, where j may take on any value up to 
t, this requires 4 x 75 x 75 x 38 = 755,000 computer words, 
beyond the core storage capacity of all but the largest com- 
puters. Since determination of the probability distribution 
of the number of overhauls is the only advantage to computing 
n(K,t,J,X), and since, as we have seen, we can compute the 
expected value of this random variable by other methods, under 
the circumstances omission of this calculation seems a reason- 
able trade-off, provided we assume that decisions will be risk 
neutral (i.e., based on expected cost). Given a computer with 
sufficient capacity, calculation of m(k,t,j,x) would be no 
more difficult conceptually than it was in the earlier model. 

The state probabilities (P(k,t,j)) resulting from the 
model are calculated utilizing the transition probabilities of 
Table 4 for a demand dependent overhaul policy. Because of 
the large number of states, these probabilites have been 


aggregated over all j for a given pair, k,t. This results in 
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the probability of the ship's being in condition k at age t 
without consideration of the time since last overhaul. This 
probability will be denoted by $(k,t). Values Of QUE t) ROT 
this example are displayed in Appendix F. 

The expected number of overhauls the ship will experience 
over its life is 5.546. This value was obtained in the same 
manner as that for the example in Chapter II: The probability 
of a ship's being in the "demand" condition (condition 4, in 
this case) at any given ship age is the probability of being 
in overhaul at that age. Since only one time unit is 
involved, this probability is also the expected number of 
overhauls that the ship will undergo at that age. The 
expected number of overhauls throughout the ship's life is 
obtained by summing these values over the ship's life, i.e., 
all possible ages. 

In the following chapter, a comparison of the demand 
dependent overhaul policy reflected here will be made with a 


policy of some fixed overhaul cycle length. 
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Е. RELAXATION OF CERTAIN ASSUMPTIONS: 


Throughout the development of the various models, we have 
retained certain assumptions: That, under a demand dependent 
policy, a ship would not be overhauled until it reached the 
worst condition; that, when it reached this condition, it 
would be overhauled with probability 1; and, that this over- 
haul would return the ship to condition 1 with certainty. 
While we will continue to make these assumptions, the approach 
can accomodate relaxing them with little difficulty. Since 
being in the worst condition would no longer be tantamount 
to being in overhaul, some additional computation would be 
required to determine the probability of being in overhaul, 
given the ship's age and the length of time since its last 
overhaul. Further, since overhaul would no longer ensure 
the ship's returning to condition 1, P(1,t,l) would no 
longer be equal to unity (for all t), and one or more of 
P(k,t,1) would have values greater than zero for k = 2,3,4. 
Again, the change would be more computational Seng concep- 


tual. 
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IV. A COMPARISON OF OVERHAUL POLICIES 
A. A SCHEDULED OVERHAUL POLICY: 


We will now use the most comprehensive model to consider 
a policy by which a ship is overhauled at some specified 
interval, rather than waiting for some demand condition to be 
met before overhauling. This is the policy currently in 
effect in the Navy. To do this requires very little modifi- 
cation to this model. For an overhaul cycle length (i.e., 
end-of-overhaul to end-of-overhaul) of j time units, we modify 


p(kli,t,j-1) as follows: 


p(k|4,t,j-1) = 1, k= 


| 
= 
~ 
~ 
~ 
ОЈ 


p(k|i,t,j-1) 


|| 
о 
s 
x 
|| 
E 
N 
> 
H 
|| 
pa 
N 
~ 
Lo 


In this manner we artificially "force" a demand condition on 
the ship. The mechanics of the model do not require any 
change at all. This modification does reflect the assumption 
that, if the ship should reach condition 4 prior to the end of 
the cycle, overhaul would take place ahead of schedule with 
the next cycle beginning at that time. 

We may use the multi-condition, age and time-since- 
overhaul dependent model, then, to compare policies of over- 
hauling at some fixed interval with the demand dependent pol- 
ciy described earlier. We will utilize our previous example 


(Section D of Chapter III) to demonstrate this. 
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With the model modified to reflect a fixed overhaul pol- 
icy as described above, the results shown in Table 5 were 
obtained. In the example of Chapter II we used the highest 
probability of being in overhaul over the ship's life as a 
basis for determining the necessary fleet size to meet speci- 
fied requirements. This highest probability is shown in Table 
5, along with the corresponding value of j, as well as the 
cycle length (in years) and the expected number of overhauls 
the ship will experience under that policy. The symbol "e" 
implies a strictly demand dependent policy. We observe that 
the highest probability of being in overhaul occurs at the 
end of the first cycle in each case (i.e., at j = 3 for a one 
year cycle length, j = 6 for a two year cycle length, etc.). 
To use these values to determine a required fleet size would 
imply a whole fleet of exactly the same age. Although we 
accepted this assumption in the example of Chapter II, such 
an assumption is obviously unrealistic and would distort the 
decision making process in this case. Even if all the ships 
in a fleet were new at the same time, we must realize that in 
actuality the decision maker would stagger the point at which 
ships' overhaul cycles were begun. Consequently, we must find 


Some way to overcome this problem. 
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Expected Highest Ship Age at 

Overhaul Number of Probability which Highest 

Cycle Length Overhauls of Being in Probability 
(Years) (Lifetime) Overhaul Occurs 
1 25267 .879 3 
2 13.96 .811 6 
3 10.24 E765 9 
4 8.43 . 641 12 
5 7.42 .480 15 
6 6.82 .371 18 
7 6.42 .294 2 
8 6.16 238 24 
9 5.98 .199 27 
10 5.85 #1 71 30 
со 5.55 -121 2 

where j is the ship's age in terms of the basic time unit. 


TABLE 5 


Expected Number of Overhauls and Highest 
Probabilities of Being in Overhaul for 
Various Overhaul Policies 
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E IHE MIXEDFLEET: 


It is much more likely that in reality only a few new 
ships would be acquired at a time. We will assume that our 
ship acquisition policy calls for purchasing some number of 
new ships every two years. If we assume that this takes place 
over, say, ten years, then we essentially have five different 
types of ship to consider, since the probability of being in 
overhaul varies with ship age, even though we still assume 
that the ships are otherwise identical. 

In actual practice in the Navy a fleet of ships is con- 
tinuously experiencing replacement of old ships with new ones. 
Consequently we should examine our "mixed" fleet over some 
period where we have essentially reached a "steady state" in 
number of ships (rather than some point where we are either 
Still acquiring ships on net or have begun decreasing the 
total fleet size), say a five year period beginning when our 
ships' ages vary from 12 to 20 years. Ме will wish to examine 
the mixed fleet at perhaps five points in time during this 
period, say a year apart, in order again to determine the 
total fleet size required to provide us with a given number of 
operational ships with a specified probability.  Unfortu- 
nately, since all our ships do not simultaneously have the 
same probability of being in overhaul we cannot apply the 


binomial probability mass function as we did in Chapter II. 





53 

There is, however, another method of approaching the 
problem. First, assume that we have a relatively large num- 
ber of ships. Secondly, suppose that the probability of one 
ship's being in overhaul is independent of the probability of 
any other ship's being in overhaul. This should in general 
be true except for combatant situations, natural disasters, 
etc., or unless the capacity of the overhauling shipyards 
becomes a constraint. 

If these conditions are satisfied then we may apply the 
central limit theorem of probability. Essentially, the cen- 
tral limit theorem states (for our application) that the sum 
Of a sequence of N independently distributed random variables 
converges in distribution as N>w to a random variable that is 
normally distributed with mean equal to the sum of the com- 
ponent means and with variance equal to the sum of the com- 
ponent variances, provided that Liapunov's condition is met. 
Cramér* states Liapunov's condition as follows: Let 21, 
гот... be independent random variables, and denote by m. and 
m the mean and the standard deviation of zi. Suppose that 


the third absolute moment of 2. about its mean 


Cu E 3 


is finite for every i, and write 


"See Harald Cramér, Mathematical Methods of Statistics, 
Princeton University Press, Princeton, NJ, 1946, pp. 215-218. 
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Er the condition 


lim фе. 
= 0 


k >% 
k 
is satisfied, then the sum z = | 2. is asymptotically normal 
1=1 
(m,o) where m and o are given by: 


m = my + m +... + т, 


2 2 2 2 
О = + + seo + Об“, 
AO k 


For a set of Bernoulli trials 


3 — _ - — 2 2 
p: =wE(|z, - p,1%) = pay (ej + qi) 


с ра; 
where 
d; = B Pj” 
Therefore 
k 
p> < | Р 4 - py). 
i=l 
Now: 
Е k 
с p SD 
1=1 
Therefore 
В Е -1/6 
mds В. Ти... 
с 151 1 Ж 
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со 


ТЕ the series ) р; (1 - р.) is divergent, Liapunov's condition 
т 
is satisfied and thus the variable z is asymptotically normal 


са 
P.r рар 1 
i=1 + i=l? 3 


со 


The series p р; (1 E р;), where Pi is the probability of the 
ith ship's Es in overhaul, will be divergent if for all i 
р; equals neither 0 пог 1. This can be reasonably assumed for 
our application. 

The approximating normal a will then have a 
mean of „2 Pi and a variance of LB = р;). 

Since we have our ships in five groups and are assuming 
that all ships within a group are identical, we may simplify 
the summation process somewhat as follows: If we have n; 
ships with probability p; of being in overhaul PË some speci- 
fied EI 1 = 1,...,5, then we may write: u =) п.р}, апа 
E - ) п.р; (1 - pi); in order to obtain the me. and the 
E. S“, of our approximating normal distribution. 

If we take N = Ds then, using the normal distribution, 
we may approximate T probability of having k out of N ships 


in overhaul, where k is equal to or less than a specified 


value b, as follows: 


Prik < b} го је 1⁄2 = и 


ЕУ) is defined as the cumulative distribution function for 


the unit normal probability density function, and 
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| ету /2 ду. 
The +1/2 term is added to the argument of $ in order to 
reflect the discreteness of the distribution we are approxi- 
mating. The reasoning here is similar to that in the 
DeMoivre-LaPlace limit theorem. 

Using this approach, then, we may determine the size of 
fleet necessary to meet our requirement. Obviously, given 
that we have initially assumed a fleet size as before, any 
increase or decrease may be implemented among our five dif- 
ferent groups of ships in various ways. Where we decide to 
make such an increase or decrease may have some effect on our 
outcome. 

Now that we have developed the approach to the problem 
using the normal approximation, there is no reason why we must 
restrict ourselves to ships which are identical except for 
age. We are equipped to consider ships with entirely differ- 
ent sets of material condition transition probabilities. For 
each type we must, of course, go through the procedure demon- 
strated in the example of Chapter III for determining the 


state probabilities. 





57 
C. AN EXAMPLE: 


Returning to the results of the example of Chapter III: 
Suppose we wish to consider a fleet of five different groups 
of such ships, initially varying in age from twelve to twenty 
years, over a five year period. Define T as the epoch in time 
under consideration. For convenience let T be equal to the 
age of the oldest group of ships. Table 6 shows values for т 
and the ages of the different groups, in terms of the basic 
(four month) time unit, at each of our five observations. 

For the example we will assume that the decision maker 
requires a probability of .95 for the event {20 or more ships 
not in overhaul). We will initially assume a required fleet 
size of 25 ships: five in each group. If it is necessary to 
increase the fleet size to meet the requirement, we will begin 
by adding a ship to Group 5, then one to Group 4, etc. ТЕ, оп 
the other hand, we wish to decrease the fleet size because the 
assumed size is more than sufficient to meet the requirement, 
we will remove one ship from Group 1, then one from Group 2, 
БЕС. 

In the example we will consider overhaul cycle lengths of 
four, six, eight, and ten years as well as the demand depen- 
dent case. For each of these we will take the observation 
with the largest fleet size necessary to support the decision 
maker's requirement and compare the policies on a mean value 


basis, economizing total cost. Here we will assume, as we 
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Age of 
Group 1 Age of Age of Age of AJ ОЕ 
Observation 1 (oldest) Group 2 Group 3 Group 4 Group 5 


1 60 60 54 48 42 36 
2 63 63 57 51 45 39 
3 66 66 60 54 48 42 
4 69 69 63 57 51 45 
5 72 72 66 60 54 48 


Age is given in terms of the basic time unit. 


TABLE 6 


Ship Ages for the Mixed 
Fleet Example 
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have throughout the thesis, that our decision maker is not a 
risk averter and will therefore base his decision on expected 
cost. 

Table 7 shows the values of р; (the probability of being 
in overhaul for the ith group) for each of the five observa- 
tions under each of the five overhaul policies. Using these 
values we now compute the number of ships in each group (under 
the decision rules previously established) needed to meet our 
requirement. The computer model used to do this is shown as 
Appendix G. This model, in the form shown, actually computes 
the probabilities for the event {20 or more ships not in over- 
haul} for all combinations of four, five, and six ships in 
each of our five groups, so the decision maker could examine 
various rules for increasing or decreasing the initially 
assumed fleet size. 

The requirements resulting from the computation are tabu- 
lated in Table 8 for each of the five overhaul policies we have 
chosen to consider. The n, are the requirement for the number 
of ships in each of the five groups with N being total fleet 
size. 

For the various policies, the maximum N (upon which we 


wıll base our decision) is summarized as follows: 





Four Year Overhaul Cycle 


Observation PE ча P3 
1 .11683 .10314 . 12567 
2 .10860 SEKR .11084 
3 .10600 .11683 .10314 
з .10830 .10860 .10721 
5 . 11260 .10600 „11683 
Six Year Overhaul Cycle 
Observation em Si Рз 
1 . 08313 . 09742 . 08469 
2 .08511 .09094 .08631 
3 .08619 .08313 . 09742 
4 .08656 .08511 .09094 
5 . 08980 ‚08619 ‚08313 
TABLE 7 


409749 
.10418 
912509 
.11084 


.10314 


.07609 
· 08047 
‚08469 
. 08631 


.09742 


Probabilities of Being in Overhaul at Each 
Observation for the Mixed Fleet Example 
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Ps 


.17841 
DU 
. 09749 
.10418 


. 12567 


.15261 
.10574 
. 07609 
.08047 


. 08469 


Observation 


Observation 


Eight Year Overhaul Cycle 


P1 


. 07856 
.08183 
.07494 
207598 


QE 


Ten Year Overhaul 


Рі 


07078 
.07114 
„07252 
‚07466 


‚07064 


P2 


-07533 
.07668 
‚07856 
„08183 


07494 


Po 


06968 
07034 
. 07678 
.07114 


„07252 


P3 


. 08715 
. 07215 
07533 
.07668 


.07856 


Сус1е 


Рз 


. 06854 
„06909 
„06968 
. 07034 


„07678 


TABLE 7 (continued) 


.06887 
07237 
.08715 
. 017215 


. 07533 


Pa 


06591 
. 06818 
.06854 
. 06909 


.06968 
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.10810 


.10670 


. 06887 
07217 


.08715 


.10064 
.10061 
„06591 
„06818 


‚06854 





Observation 


Demand Dependent Policy 


Py 


. 06395 
.06395 
„06395 
. 06595 


-06395 


Во 


‚06395 
„06395 
206395 
206395 


. 06395 


Рз 


.06394 
206595 
206395 
06395 


- 06395 


TABLE 7 (continued) 


Pa 


. 06165 
. 06407 
‚06394 
406395 


. 06395 


62 


„09980 


‚09980 


‚06165 


‚06407 


‚06394 





6 3 


Four Year Overhaul Cycle 


Observation E Ша но 14 Si N 
т 5 5 5 5 6 26 
2 5 5 5 5 5 25 
3 5 5 5 5 5 25 
4 5 5 5 5 5 25 
5 5 5 5 5 5 25 

Six Year Overhaul Cycle 

Observation “1 ша ES 14 B5 N 
й 5 5 5 5 5 25 
2 4 5 5 5 5 24 
3 4 5 5 5 5 24 
4 4 5 5 5 5 24 
5 4 5 5 5 5 24 

TABLE 8 


Fleet Size Requirements for the 
Mixed Fleet Example 





Observation 


Observation 


Eight Year Overhaul Cycle 


ut 2 3 
4 5 5 
4 5 5 
4 5 Б 
4 5 5 
4 5 5 


1 2 3 
4 5 5 
4 5 5 
4 5 5 
4 5 5 
4 5 5 


TABLE 8 (continued) 
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је 


24 


24 


24 


24 


24 


e 


24 


24 


24 


24 


24 
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Demand Dependent Policy 


Observation = n2 E Fa 15 NO 
1 4 5 5 5 5 24 
2 4 5 5 5 5 24 
3 4 4 5 5 5 23 
4 4 4 5 5 5 28 
5 4 4 5 5 5 23 


TABLE 8 (continued) 
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POLICY “пах 
Four Year Overhaul Cycle 26 
Six Year Overhaul Cycle 25 
Eight Year Overhaul Cycle 24 
Ten Year Overhaul Cycle 24 
Demand Dependent Policy 24 


Recalling the expected number of overhauls from Table 5, 
we may compute total expected cost for the fleet, under each 


policy, as follows: 


Total expected cost = N [Ship cost + 


Overhaul cost (Expected number of overhauls)] 


As before, we will assume a ship cost of $40 million and an 
overhaul cost of $1.5 million. Thus for the four year overhaul 


cycle policy, total expected cost would be computed as follows: 


Total expected cost 26[$40 x 10°+ $1.5 x 10° (8.43)] 


S1.37 billion. 


For the remaining policies considered, total expected cost is: 


POLICY TOTAL EXPECTED COST 
Six Year Overhaul Cycle S1.26 billion 
Elght Year Overhaul Cycle $1.18 billion 
Ten Year Overhaul Cycle $1.17 billion 


Demand Dependent Policy s I1.16 pillion 
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The decision maker would then choose the demand dependent pol- 
icy as that which would be most economical (based on expected 
costs). As one would expect, and as we see, due to the fact 
that Ben cost is relatively large compared to overhaul E 
the total expected cost is more sensitive to required fleet 
size than to expected number of overhauls. 

Suppose now that the decision maker is willing to relax 
his requirement to a probability of .90 for the event {20 or 
more ships not in overhaul}. In the same manner as before, we 


obtain the following results: 


POLICY N TOTAL EXPECTED COST 
Four Year Overhaul Cycle 25 $1.32 billion 
Six Year overhaul Cycle 24 $1.21 billion 
Eight Year Overhaul Cycle 24 $1.18 billion 
Ten Year Overhaul Cycle 23 $1.12 billion 
Demand Dependent Policy '23 STCTI Da NS 


On the other hand, if he requires a probability of, say, .98 
for the event {20 ог more ships not in overhaul], we obtain 


the following: 


POLICY N TOTAL EXPECTED COST 
Four Year Overhaul Cycle 27 $1.42" billion 
Six Year Overhaul Cycle 26 Lil lion 
Eight Year Overhaul Cycle 25 опе Бран Топ 
Ten Year Overhaul Cycle 24 SI 1 con 


Demand Dependent Policy 24 $1.16 billion 





68 


We observe that, while the decision is relatively insensitive 
to such a change in requirements, the required size of fleet 
NS not. 

Up to this point we have assumed that the cost of a 
scheduled overhaul (under a policy calling for an overhaul 
cycle of some specified length) is the same as that of an 
unscheduled overhaul (i.e., one resulting from the ship's 
reaching the demand condition prior to the end of the overhaul 
cycle). For a strictly demand dependent policy, of course, 
all overhauls are unscheduled. An unscheduled overhaul should 
be expected to cost more than a scheduled one, due to uncer- 
tainty as to the number of ships in overhaul at any given 
time, the inability to plan specific work ahead of time based 
on an individual ship, etc. 

Using our multi-condition, age and time-since-overhaul 
dependent model, the expected number of both scheduled and 
unscheduled overhauls which a ship will experience, under a 


policy calling for specified cycle lengths, can be computed as 


follows: 
T-1 
E(number of scheduled overhauls) = ў P(4,t,L.); 
t-1 | 
TL Пета 
E(number of unscheduled overhauls = ) TEE t, 3) 
t=1 j=l 


where Ly is the specified overhaul cycle length. These values 


have been computed for our example and are shown in Table 9 





Overhaul 
Cycle Expected Number 
Length of Scheduled 
(Years) Overhauls 
1 21.70 
2 AO 
3 O 
4 3.48 
5 2.27 
6 1.54 
y 1.06 
8 0.74 
9 0.52 
10 0237 
oo 0.00 


TABLE 9 


Expected Number of Scheduled and 


Expected Number 
of Unscheduled 
Overhauls 


3:97 


4.01 


4.65 


4.95 


5.42 


5.46 


Unscheduled Overhauls for the 
Multi-Condition, Age and Time- 
Since-Overhaul Dependent 


Model Example 
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for overhaul cycle lengths of one to ten years and for the 
demand dependent policy. 

Suppose we assume that the cost per overhaul of an 
unscheduled overhaul is $2 million, while that of a scheduled 
overhaul is $1.5 million. We will again examine the total 
expected fleet cost, taking this assumption into account. ТЕ 
we require a probability of .90 for the event {20 or more 
ships not in overhaul} (denote this as рез ZUR we obtain 


the following results: 


FOLICY N TOTAL EXPECTED COST 
Four Year Overhaul Cycle 25 S1.38 billion 
Six Year Overhaul Cycle 24 СаО топ 
Eight Year Overhaul Cycle 24 SI.25' billion 
Ten Year Overhaul Cycle 23 S1.18+ billion 
Demand Dependent Policy 23 g Ig pillion 


For a ВЕ (20) of .95 we obtain: 


POLICY N TOTAL EXPECTED COST 
Four Year Overhaul Cycle 26 $1.43 billion 
Six Year Overhaul Cycle 25 21. За Lion 
Eight Year Overhaul Cycle 24 ор: Zoi о 
Ten Year Overhaul Cycle 24 51.24 billion 
Demand Dependent Policy 24 s1.23 pillion 


And, for ER - „98: 
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РОГТСУ N TOTAL EXPECTED COST 
Four Year Overhaul Cycle 27 S1.49 billion 
Six Year Overhaul Cycle 26 s1l.37 pillion 
Elght Year Overhaul Cycle 25 S1.30 billion 
Ten Year Overhaul Cycle 24 S1.24 billion 
Demand Dependent Policy 24 sL.23 Тол 


As we can see, the decision as to the most economical policy 
would remain unchanged with the assumption that an unscheduled 
overhaul costs $2.0 million ($0.5 million more than a sched- 
uled one). 

We may approach the trade-off between scheduled and 
unscheduled overhauls from another viewpoint, however. Sup- 


pose we wish to find some coefficient c, where 


_ cost of unscheduled overhaul 
DE о a o eae 


cost of scheduled overhaul 


such that the decision maker would be indifferent between two 
alternative overhaul policies. To illustrate the manner in 
which we can do this, consider the ten year overhaul cycle 
policy and the demand dependent policy for Bel) = „95. 
Since the required fleet size is the same, we need only con- 
sider single ship expected overhaul cost in order to determine 


a value for c. We do this as follows: 


E(number of scheduled overhauls for 10 yr. policy) 


+ c x E(number of unscheduled overhauls for 10 yr. policy) 
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= E(number of scheduled overhauls for demand dependent 
policy) + c x E(number of unscheduled overhauls for 


demand dependent policy). 
Substituting the values from Table 9 we get: 


=> + 5.486 — 21004175. 956 


257 


с = 7g; ^ 3 


This implies that for the particular set of transition proba- 
bilities assumed for the example, the cost of an unscheduled 
overhaul would have to be of the order of five times that of 
a Scheduled one in order for the decision maker to change his 
decision to adopt the demand dependent overhaul policy. 

In our example the demand dependent overhaul policy has 
consistently been the choice, economically. This will always 
be true unless we impute a penalty to the cost of an unsched- 
uled overhaul. Given such a penalty, and different behavior 
of transition probabilities from our example (especially if 
probability of "failure" continues to increase as ship age 
and time since overhaul increase), we could easily have 
reached a different decision. 

We have shown, then, how our model can be used to compare 
alternative overhaul policies by considering total expected 
cost of a fleet of ships in addition to the cost differential 


between scheduled and unscheduled overhauls. 
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V. RECOMMENDATIONS, SUMMARY, AND CONCLUSIONS 
A. RECOMMENDATIONS FOR FURTHER WORK: 


The next step, and the most important toward actual 
implementation of the models, and in particular the most com- 
prehensive (multi-condition, age and time-since-overhaul 
dependent) model, is determination of reasonable estimators 
for transition probabilities based on the best data available. 
Once these have been determined, the model implemented, and 
transition probabilities refined so that the results of the 
model are relatively reliable, the model may be used to exa- 
mine a large number of possible maintenance strategies. For 
example, the decision maker might want to consider the alter- 
native of assigning his ships to shipyards upon their reaching 
condition 3 for limited maintenance not actually amounting to 
overhaul (as we have used the term here). There are no doubt 
other maintenance policies of this nature that can be eval- 
uated using this multi-condition, time-since-overhaul дереп- 
dent model as a basis. 

ТЕ it were found to be reasonable to describe transition 
probabilities analytically in terms of ship age and time since 
last overhaul, the model certainly could be modified to accom- 
modate such a change. As a matter of fact, transition proba- 
bilities could be computed internally by the computer model. 


This would eliminate the necessity of reading in (in some 





74 
manner) long lists of transition matrices, as we have done 
here. 

Efforts to develop criteria for an acceptable material 
condition index for ships should continue, although the com- 
plexity and variety of configurations encountered in warships 
does not make this an easy assignment. An index of some kind 
is essential if one is to expect valid results from this 
model. Іп the interim, however, the model should be imple- 
mented utilizing whatever rough index of material condition 
that can be developed in the short run. 

It must be emphasized that results obtained from the 
models are only as good as their input. The models are geared 
to categorization of material condition either by subjective 
judgment or by some quantitative scheme for determining a 
material condition index. Critical to the model's usefulness 
as an aid to decision making is accurate determination of 
transition probabilities, as indicated above. As pointed out 
in Appendix D, these should be estimated initially from the 
best historical data available, using criteria which, once 
decided upon, should be applied consistently to all ships. 
Then, using the Bayesian approach,? transition probabilities 
should be updated as actual ship behavior is observed. 

“For one example of hov this approach can be used, see 
Chapter 5 of John W. Devanney III, Marine Decisions Under 
Uncertainty, National Sea Grant Program, Sea Grant Project 


GH-88, M.I.T., Cornell Maritime Press, Cambridge, Maryland, 
November 1971. 
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We have demonstrated the model using only one Set of 
transition probability matrices. It would be useful to hypo- 
thesize transition probabilities based on different assumed 
forms of behavior of the ship's material condition (e.g., 
deterioration more or less linear over time, etc.) and compare 
the results with those obtained here. In this manner some 
"feel" may be obtained for the dynamics of the model. 

An approach to modeling the change in a ship's condition 
over time that deserves some consideration involves the use of 
a semi-Markov process model. In a semi-Markov process, the 
states occupied on successive transitions are governed by the 
transition probabilities (as in our model) of a Markov process 
(called the imbedded Markov process). The time the process 
stays in a particular state, however, is a random variable 
(either discrete or continuous) described by an exponential 
density (Poisson process). This random variable is referred 
to as the "holding" time and is dependent on the state pres- 
ently occupied and on the state to which the next transition 
will be made. Howard, Volume II,° provides a rather thorough 


discussion of semi-Markov processes. 


SRonald A. Howard, Dynamic Probabilistic Systems, Volume 


en — —T s sp e 


New York, 1971, Chapters 10 and 11. 
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B. SUMMARY AND CONCLUSIONS: 


Beginning with a very simple description of the behavior 
of a ship's material condition, we have developed a set of 
models which, while still conceptually simple, is capable of 
handling a rather wide range of ship deterioration behavior. 
The approach used in developing the models considers the ship 
as a single component, rather than a set of many smaller com- 
ponents. Modeling the behavior of the ship's material con- 
dition over time by a transient Markov process represents an 
unusual if not unique approach to this area of no small con- 
cern to the Navy. 

With the central limit theorem, as we have seen, there is 
no problem in going from considerations of a fleet of identi- 
cal ships, all of the same age, to a "mixed" fleet, where not 
only are ships not all of the same age but are dissimilar in 
other respects as well. 

We have seen how the models can be used to determine the 
effect of various overhaul policies on ship acquisition 
requirements as well as on total expected fleet cost. Fur- 
ther, we have seen the importance of considering the real cost 
differential between scheduled and unscheduled overhauls. 

The significant conclusions to be drawn from this thesis 


are threefold: 


(1) That the use of transient Markov processes to model 


ship deterioration and failure behavior is reasonable. 
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(2) That such models are potentially of significant use 


to Navy (and other) decision makers. 


(3) That a great deal of initial effort is required in 
the area of data analysis to develop material condition 
indices and reasonably accurate estimates of transition prob- 


abilities. 
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0. 500000 
0.100000 
0.910000 
0. C90000 
0.909000 
0.091000 
0. 509100 
0.090900 
0. 509090 
0.090910 
0.909091 
0. С9 0909 
0.818181 
0.181818 
0. #36363 
0.163636 
О. 832 726 
0.167273 
0.833454 
0. 166545 
0. 833308 
0.166691 
0. 833337 
0. 1666€2 
0.749997 
0. 250001 
0. 774999 
0.224999 
0. 767499 
0.232500 
0.769749 
0. 230250 
0.769074 
0. 230925 
0.769276 
0.230722 
О. 692288 
0.307710 
0.723083 
0. 276915 
0.710765 
0.289233 
0. 715652 
0.284306 
0.713721 
0. 2862 77 
0.714509 
0.285488 
О. 642743 
О. 357255 
0. 678626 
0.321371 
0.660684 
0. 339313 
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0.669655 
0. 330342 
0.665170 
0.334828 
0. 667412 
0.332585 
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©` С^ с^ © с^ m. QO. OL ORO OO RO ЛА АЛЛА АЛАЛАЛ АЛ IR Dia Vin vn + + + ES PS P P bb va Ww ww Ыыы ыш ммм м 


Мо ме N ке IN) e N = N) = N = IN) m N = N = Кр ре Кр ра Кр ме (у ке ND) ке IN) m N ra IN) m IN); ке ру ке Ку МО Fa y N = N = N A N = N = lo^ 


Моне нели ПОВИ о АИ ОЈО Ооо дого пукнала го каа e na ra ar T 


0.900000 
1.000000 
0.100000 
с. 000000 
Ce 910000 
0.900000 
С. 900000 
1.000000 
0.090000 
0. 100000 
0.100000 
C. 000000 
C.909000 
0.910000 
C. 910000 
0.900000 
0.900000 
1. 000000 
0.091000 
С. 090000 
С. О90000 
0.100000 
C.100000 
0.000000 
0.909100 
C. 909000 
0.909000 
Ce 910000 
Ce 310000 
0.900000 
C. 900000 
1.000000 
0.090900 
0.091000 
0.091000 
С. 090000 
0.090000 
0.100000 
С. 100000 
0.000000 
0.909090 
С. 909100 
0.909100 
С. 909000 
Ce909000 
0.910000 
Ce 910000 
0.900000 
0.900000 
1. 000000 
0.090910 
C. 090900 
C. 090900 
C.091000 
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C.091000 
C. 090000 
0.090000 
Ce 100000 
0.100000 
0.000000 
С. 818182 
0.818180 
0.818180 
0.818200 
С. 818200 
С. 818000 
0.818000 
0.820000 
0.820000 
0.800000 
Ce 800000 
1.000000 
0.181818 
0.181820 
0.181820 
0.181800 
С. 181800 
0.182000 
0.182000 
0.180000 
0. 180000 
С. 200000 
0.200000 
Ce 000000 
Ce 836363 
0,836364 
0.836364 
0.836360 
0.836360 
0. 836400 
0.836400 
0.836000 
0.836000 
0.840000 
C. 840000 
0.800000 
0.800000 
1. 000000 
0.163636 
C.163636 
0.163636 
0.163640 
0.163640 
0.163600 
0.163600 
0. 164000 
0.164000 
0.160000 
0.160000 
0.200000 
0. 200000 
0.000000 
0.832727 
С. 832727 
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0.832727 
0.832728 
0.832728 
0.832720 
0.832720 
0.832800 
0. 832800 
0.832000 
0.832000 
C. 840000 
0.840000 
C. 800000 
0.800000 
1.000000 
с. 167273 
0.167273 
0.167273 
С. 167272 
0.167272 
0.167280 
0.167280 
0.167200 
C. 167200 
0.168000 
0.168000 
0. 160000 
0.160000 
6.200000 
0.200000 
0.000000 
Ce 833454 
0.833454 
0.833454 
С.833454 
0.833454 
С. 833456 
0.833456 
0.833440 
0.833440 
0.833600 
0.833600 
0.832000 
0.832000 
0.840000 
C.840000 
0.800000 
C. 800000 
1.000000 
0.166545 
С. 166545 
0.166545 
С.166545 
0.166545 
0.166544 
С. 166544 
0.166560 
0.166560 
0. 166400 
0.166400 
6.168000 
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0.168000 
Ce 160000 
0.160000 
0.200000 
С. 200000 
0.000000 
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0.900000 
0.000000 
0.000000 
0.100000 
0.810000 
0.000000 
0.1C0000 
0.090000 
0.000000 
О. 060000 
0.729000 
0. 060000 
0.180000 
0.081000 
0. 060000 
0.010000 
0.000000 
0.000000 
0.656100 
0.000000 
0.243000 
0.072900 
0.010000 
0.018000 
0.000000 
0.0C0000 
0.000000 
0.000000 
0. 550490 
0.000000 
0.291600 
0. 065610 
0.027000 
0.024300 
0.0C0000 
0.001000 
0.000000 
0.000000 
0.000000 
0.000000 
0.531441 
0.000000 
0.328050 
0.059049 
0.048600 
0.029160 
0.001000 
0.002700 
0.0C0000 
0. 0C0000 
0.000000 
0. 000000 
0.000000 
0.000000 
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0.425153 
0.000000 
0.321489 
0.106288 
0.068040 
0.065610 
0.003500 
0,009720 
0.0C0000 
0.000200 
0.000000 
0.000000 
0.0C0000 
0.000000 
0.0C0000 
0.000000 
0.340122 
0.0C0000 
0.363479 
0.085030 
0.120042 
0.064298 
0.012520 
0.013608 
0.000200 
0.0C0700 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.0C0000 
0.000000 
0.000000 
0.272098 
0.000000 
0.375814 
0.068024 
0.160331 
0.072696 
0.023624 
0.024008 
0.000860 
0.002504 
0.0C0000 
0.000040 
0.000000 
0.000000 
0.000000 
0.000000 
0.0C0000 
0.000000 
0.000000 
0.000000 
0.217678 
0.000000 
0.3€8675 
0.054420 
0.2C0961 
0.075163 
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0.042907 
0.032066 
0.003192 
0.004725 
0.000040 
0.000172 
0.000000 
0.000000 
0.000000 
0.000000 
0.0C0000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.174142 
0.000000 
0. 349 36C 
0.043536 
0.235931 
0.073735 
0.066392 
0.040192 
0.007278 
0.008581 
0.000204 
0.000638 
0.0C0000 
0.000008 
0.000000 
0.0C0000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0. 000000 
0.139314 
0.000000 
0.323023 
0.034828 
0.262480 
0.069872 
0.093306 
0.047186 
0.014404 
0.013278 
0.000802 
0.001456 
0.000008 
0.000041 
0.000000 
0.0C0000 
0.000000 
0.000000 
0.000000 
0.000000 
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0. 000000 
0.000000 
0. 000000 
0.00000C 
0.000000 
0.060000 
0.097520 
0. 060000 
0.260945 
0.041794 
0. 2 53608 
0.096907 
0. 112500 
0.078744 
0.023361 
0. 027992 
0.002017 
0.004321 
0.000046 
0.000240 
0. 0C0000 
0.000002 
0.000000 
0.0C0000 
0.000000 
0.000000 
0.000000 
0.000000 
0. 060000 
0.060000 
0.000000 
0. 000000 
0.000000 
0.000000 
0.0€8264 
0.000000 
0.224455 
0.029256 
0.274432 
0. 078283 
0.157494 
0.076082 
0.044345 
0.033750 
0. 065733 
0.007008 
0.000273 
0. 000605 
0.000002 
0.000014 
0.0C0000 
0.000000 
0.0C0000 
0.0C0000 
0.000000 
0.0C0000 
0.000000 
0.000000 
0.000000 
0.000000 
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0.000000 
0.000000 
0. 000000 
0.000000 
0.047785 
0. 000000 
0.186375 
0.020479 
0.270386 
0.067337 
0.186328 
0.082330 
0.064791 
0.047248 
0.011022 
0.013303 
0.000796 
0.001720 
0.000016 
0.000082 
0.000000 
0. 000001 
0.000000 
0.000000 
0. 060000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0. 060000 
0.000000 
0.000000 
0. 060000 
0.033445 
0.0C0000 
0.150941 
0.014335 
0.256607 
0.055912 
0.212759 
0.081116 
0.092602 
0.055898 
0.021018 
0.019437 
0.002277 
0.003306 
0.000093 
0.000239 
0.000001 
0.000005 
0.000000 
0.0C0000 
0.000000 
0.000000 
0.0C0000 
0.000000 
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0. 000000 
0.000000 
0.000000 
0.0C0000 
0.000000 
0.0C0000 
0.000000 
0.000000 
0.000000 
0.000000 
0. 02341 4 
0.000000 
0.119994 
0. 010035 
0.235537 
0.045282 
0.230047 
0.076982 
0.120720 
0.063828 
0.034150 
0.027781 
0.004900 
0.006 306 
0. 000304 
0.000683 
0.000005 
0.000028 
0.000000 
0.0C0000 
0.000000 
0.000000 
0.0C0000 
0.000000 
0.000000 
0.000000 
0.000000 
0.0C0000 
0.000000 
0.000000 
0. 060000 
0.060000 
0.000000 
0. 000000 
0.000000 
0.000000 
0.016390 
0.000000 
0.094031 
0.0C7024 
0.210158 
0. 035998 
0.238015 
0.070661 
0.148332 
0.069014 
0.051686 
0.026216 
0.009736 
0.010245 
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0.000896 
0.0C1470 
0.000031 
0. 060091 
0.000000 
0.000002 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.0C0000 
0.000000 
0.000000 
0.0C0000 
0.000000 
0.000000 
0.000000 
0.000000 
0.0C0000 
0.000000 
0.000000 
0.009834 
0.000000 
0.063443 
0,066556 
0.162093 
0.037612 
0.213470 
0.084063 
0.158013 
0.095206 
0.067227 
0.059333 
0.016087 
0.020674 
0.002008 
0.003894 
0.000110 
0.000358 
0.000002 
0.000013 
0.000000 
0.000000 
0.0C0000 
0.000000 
0.0C0000 
0.0C0000 
0.000000 
0.0C0000 
0.000000 
0.000000 
0.0C0000 
0.000000 
0.0C0000 
0.000000 
0.000000 
0.0C0000 
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0.000000 
0. 000000 
0.000000 
0.000000 
0. ОС5900 
0.000000 
0.044622 
0.003934 
0.134868 
0.025377 
0.212145 
0.064837 
0.190014 
0.085388 
0.059669 
0.063205 
0.0 30326 
0.026891 
0.005099 
0.006435 
0.000424 
0.000803 
0.000014 
0.000044 
0.000000 
De 060001 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.0C0000 
0.000000 
0.0C0000 
0.000000 
0.000000 
0.000000 
0.000000 
0.0C0000 
0.000000 
0.000000 
0.0C0000 
0.000000 
0.000000 
0.000000 
0.000000 
0.003540 
0.000000 
0.030707 
0.0C2360 
0.106298 
0,017849 
0.192124 
0.053947 
0.159396 
0.084858 
0.123007 
0.076005 
0.045087 
0.039868 
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0.009494 
0.012131 
0.001058 
0.002040 
0.000052 
0.000170 
0.000001 
0.000005 
0.000000 
0. 0C0000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0. 000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0. 000000 
0.000000 
0.000000 
0. 000000 
0.000000 
0.002124 
0.0C0000 
0.020784 
0.001416 
0.081627 
0.012283 
0.169222 
0 „042515 
0.204496 
0.076850 
0.149809 
0.079758 
0.066920 
0.049203 
0.017827 
0.018035 
0.002674 
0.003798 
0.000201 
0.000423 
0.000006 
0.000021 
0.0C0000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
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0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0. 000000 
0.000000 
0.000000 
0. 0C0000 
0.000000 
0.000000 
0.001274 
0.000000 
0.013887 
0.000850 
0.061259 
0.008314 
0.144052 
0.032651 
0.199547 
0.067689 
0.169644 
0.081798 
0.089354 
0.059924 
0.028731 
0.026768 
0.065402 
0.007131 
0. 060544 
0.001070 
0.000024 
0.0C0080 
0.000000 
0.000002 
0.000000 
0.000000 
0.000000 


0.000000 


0.000000 
0.0C0000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.060000 
0.000000 
0.000000 
0.000000 
0.000000 
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0. 060000 
0.000000 
0.000000 
0.060000 
0.000765 
0. 060000 
0.009182 
0.000510 
0045065 
0.005555 
0.119082 
0.024504 
0.187417 
0.057621 
0.183585 
0.079819 
0.113536 
0.067858 
0.044006 
0.035742 
0.010372 
0.011492 
0.001396 
0.002161 
0.000095 
0.000217 
0.000003 
0.000010 
0.000000 
0.0C0000 
0.000000 
0.000000 
0.0C0000 
0.000000 
0.000000 
0.000000 
0.000000 
0.060000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0. 0C0000 
0.000000 
0.000000 
0.0C0000 
0.000000 
0. 060000 
0.000382 
0.000000 
0. 005101 
0.000382 
0.028089 
0.004591 
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0.084045 
0.022535 
0.151329 
0.059541 
0.171611 
0.093708 
0.124626 
0.091792 
0.057745 
0.056768 
0.016678 
0.022003 
0.002859 
0.005186 
0.0C0265 
0.000698 
0.000011 
0.000048 
0.000000 
0.000001 
0.000000 
0.000000 
0.0C0000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.0C0000 
0.000000 
0.000000 
0.0C0000 
0.000000 
0.000000 
0.000000 
0.000000 
0.0C0000 
0.000000 
0.000000 
0.0C0191 
0.000000 
0.002933 
0.000191 
0.018635 
0.002550 
0.0 64557 
0.014045 
0.135206 
0.042022 
0.179514 
0.075665 
0.154105 
0. 085806 
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0.085641 
0.062313 
0.030342 
0.028872 
0.006615 
0.0C8339 
0.0C0830 
0.001429 
0.000053 
0.000132 
0.000001 
0.000006 
0.000000 
0.000000 
0.000000 
0.000000 
0.0C0000 
0.000000 
0.000000 
0. 0C0000 
0,000000 
0.000000 
0.000000 
0.000000 
0.0C0000 
0.009000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.0C0000 
0.000000 
0.000000 
0.000000 
0.000000 
0.0C0000 
0.000000 
0.000000 
0, 0C0096 
0.000000 
0.001657 
0.000096 
0.011868 
0.001486 
0.046323 
0.009318 
0.109625 
0.032278 
0.165422 
0.067603 
0162858 
0. 089757 
0.165133 
0.077053 
0.044044 
0.042820 
0.011647 
0.015171 
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0.001845 
0.003308 
0.000159 
0.000415 
0.000006 
0.000027 
0.0C0000 
0.000001 
0.000000 
0.0C0000 
0.000000 
0.000000 
0.000000 
0.000000 
0.0C0000 
0.000000 
0.000000 
0.000000 
0.000000 
0.0C0000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.0C0000 
0.000000 
0.0C0000 
0.000000 
0.000000 
0.0C0000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000048 
0.000000 
0.0C0924 
0.000048 
0.007400 
0.0C0829 
0,0 32475 
О. 065934 
0.087091 
0.023161 
0.150314 
0.054813 
0.171186 
0.082711 
0.129619 
0.081429 
0.0€4842 
0.052567 
0.020995 
0.022022 
0.004230 
0.005823 
0.000495 
0.000922 
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0.000030 
0.000080 
0.000001 
0. 0C0003 
0.000000 
0. 000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.0C0000 
0.000000 
0.000000 
0.0C0000 
0.000000 
0.0C0000 
0.000000 
0.000000 
0.0C0000 
0.000000 
0.0C0000 
0.0C0000 
0.000000 
0. 0C0000 
0.0C0000 
0.000000 
0.0C0000 
0.000000 
0.0000 24 
0. 000000 
0.000510 
0. 0C0024 
0.004529 
0.000462 
0.022174 
0.003700 
0.066707 
0.016240 
0.129969 
0.043545 
0.1£8304 
0.075157 
0, 1 46235 
0.085593 
0. 0 24988 
0.064810 
0.032519 
0.032421 
0.007938 
0.010497 
0.0C117C 
0.002115 
0.000094 
0.000247 
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0.000003 
0.000015 
0.000000 
0.000000 
0.000000 
0.0C0000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.0C0000 
0.000000 
0.000000 
0.000000 
0.000000 
0. 0C0000 
0.000000 
0.000000 
0.0C0000 
0.000000 
0.0C0000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0. 0C0000 
0.000000 
0.000000 
0.000000 
0.000000 
0. 000000 
0.000000 
0.000012 
0.000012 
0.000279 
0.000255 
0.002727 
0.002264 
0.014787 
0.011087 
0.049593 
0.033353 
0.108530 
0.064985 
0.159309 
0.084152 
0.158712 
0.073119 
0.107303 
0.042494 
0.048681 
0.016260 
0.014467 
0.003969 
0. 062700 
0.000585 
0.000295 
0.000047 
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0.000017 
0. 060002 
0.000000 
0.000000 
0.000000 
0.000000 
0. 0C0000 
0.000000 
0.000000 
0.000000 
0.000000 
0.0C0000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0. 000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.0C0000 
0,060000 
0,000000 
0. 0C0000 
0,000000 
0.000000 
0.000000 
0.000000 
0.0C0000 
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The expected number of overhauls 
the ship will experience is: 


6.37877 
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APPENDIX B 


Computer Program for Determining 
Required Fleet Size, 


Homogeneous Fleet 
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APPENDIX C 


Recursion Relationships for Age and 
Time-Since-Overhaul 


Dependent Models 
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APPENDIX C 
RECURSION RELATIONSHIPS FOR AGE AND 
TIME-SINCE-OVERHAUL DEPENDENT 


MODELS 


The following are the recursion formulas developed for 
the age and time-since-overhaul models for both 2 and 4 con- 


ditions: 
l. The 2 Condition Model: 


Plk,t,j|kort,yr3o) 9 P(1,t-1,3-1|k , t; 3 ) p (k| 13, t-1, 3-1) 
j = 2,...,t 


K 


172 


Since being in condition 2 is tantamount to undergoing an 
overhaul which will return the ship to condition 1 with cer- 
tainty, and since j takes a value of 1 only upon the ship's 


coming out of overhaul, 


-1 
) P(2,t-1,4|k9,t, јо) j=1, k=1 
=1 


t 
P(k,t,j[k rt 13,) = 2 


Also, 
T(k,t,j,x) = m(1,t-1,3j-1,x)p(k/1,t-1,3-1) 


J=2 ани ЦЕН к=] 
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T(k,t,j,x) = T(1,t-1,j-1)p(k|1,t-1,3j-1), 


J=2,...,t; k=2. 


Again, 
t 
T(k,t,),x) = 42, 


0 


n(2,t-1,2,x) j=l, k=l 


j=1, k=2. 


2. The multi-condition (4) model, following the same approach: 


3 
E Pt iket ЧЕ 


j=2,.. 


Е 


ENS Sk ,t ,j )=4: 
© -0O O gei 


0 


3 


) T (i,t-1,j-1,x)p(k]i,t-1,j-1) J 


i-1 


k=1,...,4 


}: P(4,t-1,2|k_,t 13 ) j=l, К=1 


О 


j=l; k=2,3,4. 


Sh I S 
kz1,2,3 


3 
) m(1,ta1, 71 AS E E. 


i=l 
Bc t,],x) = 
t-l 


) m(4,t-1,2,x) 


L=1 


0 


k=4 


j=1, k=1 


j=l; k=2,3,4. 
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APPENDIX D 


Determination of Transition Probabilities 
for the Example for the Multi- 
Condition, Age and Time- 
Since-Overhaul Dependent 


Model 
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APPENDIX D 
DETERMINATION OF TRANSITION PROBABILITIES 
FOR THE MULTI-CONDITION, AGE AND TIME- 


SINCE-OVERHAUL DEPENDENT MODEL 


The available data provided the dates of overhauls and 
restricted availabilities for CRUDESLANT ships from 1963 to 
1971 (approximately). Thirty-three ships were selected for 
which such data was available for a complete cycle (i.e., 
overhaul to overhaul). These ships were grouped into the 
three age categories described in Chapter III. For each cate- 
gory a histogram (Figure D-1) was developed showing the total 
number of restricted availabilities occurring during each time 
unit after overhaul (up to three years for category I and four 
years for categories II and III ~ since COMCRUDESLANT has 
never operated with an overhaul policy with cycles in excess 
of four years, little data was available for ships out of 
overhaul for longer than that period). 

The results of this were then examined for some logical 
grouping of ships based on time since last overhaul. It 
appears that in each case the number of RAV's during the first 
unit of time is markedly higher than that for subsequent time 
periods. Further, there also appears to be a relatively well 
defined difference between the second and third periods and 
later periods. Therefore, the first category of time out of 


overhaul consists simply of j = l1, the second category only of 





Total Number of RAV's 


о нм ооло мс 


O PN URB 


Or NU ANOK 


1489 


Category I 
Ship Age 1-3 Years, Inclusive 
11 Ships 





Category II 
Ship Age 4-13 Years, Inclusive 
8 Ships 


Category III 
Ship Age 14-25 Years,Incl. 


| | | 14 Ships 
1 2 3 4 5 6 7 8 9 10 11 12 


Time Units Since Overhaul 


FIGURE D-1 


Histogram of RAV Frequency 
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ПК= 2,3, while the third therefore consists G£ j = ME 
(unless j is otherwise restricted by some overhaul policy). 

With time since last overhaul categorized thusly, some 
estimate of state probabilities was required. This was done 
by first making the arbitrary assumption that one out of 
every three RAV's (which, in general, represent unexpected 
equipment failures of a relatively critical nature) was due 
to a ship's being in condition 4. Therefore, one-third the 
number of RAV's averaged over time and number of ships (for 
each category - there now being 3 x 3, or 9, categories), 
provided an estimator of the probability of being in condition 
4. It was then arbitrarily assumed that the state probabili- 
ties for condition 3 were equal to those for condition 4, and 
that those for condition 1 equaled those for condition 2. 
Since these must sum to unity, these assumptions define the 
state probabilities. Working "backwards" from the estimates 
of state probabilities, estimators for transition probabili- 
ties were developed, with the primary assumption being that 
a ship would be more likely to remain in the same condition 
over a transition than it would be to change to any other 
given condition. The matrices of these transition probability 
estimators comprise Table (4) in Chapter III. 

The high "failure rate" immediately following overhaul 
is somewhat counterintuitive. It is felt by the author that 
this is due primarily to two factors. First, when a ship is 


overhauled, disruption to major mechanical or electrical/ 
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electronic systems is frequently involved. Such disruptions 
may in some cases result in casualties to these systems once 
operation is resumed. Secondly, what is administratively 
described as a restricted availability, immediately following 
overhaul, may in fact be a continuation of the overhaul to 
complete one or more unfinished jobs. This could result from 
official disfavor toward failing to complete an overhaul on 
schedule. 

Certain of the assumptions made in developing the transi- 
tion probability matrices would not be tenable in actually 
applying the model. However, for the results of the model to 
be reasonably valid, it is important that the initial values 
estimated for transition probabilities be based on complete, 
accurate data insofar as possible, and the criteria for esti- 
mating a ship's material condition be applied in a consistent 
manner to all ships. As noted in Chapter II, while efforts 
continue to develop a material condition index which will 
permit a quantitative description of a ship's material condi- 
tion, no such usable index currently exists. Some reasonable 
approximation of such an index is required, then. Of signifi- 
cant value in developing such an approximation might be the 
results of type commanders' material inspections, or inspec- 
tions by the Board of Inspection and Survey (INSURV Board). 
Other possible parameters that could be used in estimating a 
ship's material condition include RAV frequency (as in the 


example), CASREPT data (frequency, severity, etc.), and the 
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number of man-hours expended in maintenance or the volume of 
deferred maintenance (under the "3-M" system) reported by a 
ship. Here one should be careful, however, as inconsistencies 
occur from ship to ship in the accuracy and thoroughness of 
the reporting of maintenance data. 

Whatever combination of parameters is decided upon, the 
best approach to estimating transition probabilites is prob- 
ably to choose several representative ships and use past data 
to "track" them over time, observing how they changed in 
material condition and noting age and time since overhaul at 
the time changes took place. 

It should be obvious that once transition probability 
estimators have been determined and a model such as this 
implemented, the estimators (as well as the criteria for 
grouping ships into categories) should be updated as behavior 


of ship material condition over time is actually observed. 


123 


APPENDIX E 


Computer Program for the Multi- 
Condition, Age and Time-Since 


Overhaul Dependent Model 
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APPENDIX F 


Aggregated State Probabilities (9(k,t)) 
for the Multi-Condition, Age and 
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APPENDIX G 


Computer Program for Determining 


Required Fleet Size, Mixed Fleet 
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